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Abstract 

We analyse the effect of general brane kinetic terms for bulk scalars, 
fermions and gauge bosons in theories with extra dimensions, with and 
without supersymmetry. We find in particular a singular behaviour 
when these terms contain derivatives orthogonal to the brane. This is 
brought about by 6(0) divergences arising at second and higher order 
in perturbation theory. We argue that this behaviour can be smoothed 
down by classical renormalization. 
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1 Introduction 



Field theoretical models with extra dimensions typically include one or more 
lower dimensional defects on which some fields can be localized. These de- 
fects ( "branes" from now on) can have different microscopic realizations: D- 
branes, domain walls, cosmic strings, orbifold fixed points, .... In general, 
the fields propagating in all the dimensions (bulk fields) couple to fields prop- 
agating on the branes [H |2] . Moreover, it is possible that the action includes 
terms localized on the branes involving only bulk fields ("brane terms", for 
short). They may have important consequences. A crucial observation is 
that brane terms are induced by radiative corrections. This is allowed as 
Poincare invariance is broken by the branes. Quantum divergences localized 
on the branes can appear in the correlators of bulk fields when the bulk fields 
are coupled to fields living on them Furthermore, in orbifolds this can 
occur even in the absence of brane couplings at tree level 0j. In any case, it 
is clear that the presence of divergent radiative corrections localized on the 
branes requires the introduction of divergent brane counterterms. Therefore, 
the corresponding brane couplings run with the renormalization scale, and 
cannot be set to zero at all scales. Furthermore, this indicates that they 
should be considered free parameters of the theory and be included at tree 
level (see |5] for a detailed argument). Of course, since interacting theories 
in higher dimensions are nonrenormalizable, these considerations are only 
meaningful in the framework of effective field theory, and it is important 
to take into account possible suppressions of the brane couplings by inverse 
powers of the cutoff scale A. 

Brane terms are always suppressed by 1/A n , n > 1, with respect to the 
analogous bulk terms. Therefore one would expect that their effect on ob- 
servables should be a small perturbation of the corresponding results without 
brane terms. In this article we show that, contrary to this expectation, cer- 
tain brane operators modify severely some observables — like Kaluza-Klein 
(KK) masses — even when their coefficients are very small. This breakdown 
of perturbation theory has to do with the interplay between the thin brane 
limit and the limit of large A, in which the coefficients of the brane operators 
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approach zero. Unless otherwise stated, we will always consider infinitely 
thin, delta-function branes, although we will often make use of thick branes 
as intermediate regulators. On the other hand, the underlying fundamental 
theory is usually assumed to be smooth in its parameters, so that all its low 
energy predictions should be perturbatively calculable using some effective 
theory. From this point of view, the relevant question is how to incorporate 
in a low energy effective description a mechanism to make the expansion in 
A well defined. Here we sketch a possible answer in the spirit of perturbative 
renormalization theory. 

We shall concentrate on brane kinetic terms. Kinetic terms localized on 
a brane were first invoked by Dvali, Gabadadze and Porrati j3] to recover 4D 
gravitational behaviour at short distances in infinite extra dimensions. The 
same mechanism was subsequently applied in jB] to gauge theories. In jjj, 
the one loop radiative corrections to brane terms and their first order im- 
pact on KK spectra were calculated in different theories. In j8], a detailed 
computation of infinite and finite radiative corrections to bulk and brane 
kinetic terms of gauge fields (with and without Hosotani mechanism) was 
performed. The phenomenology of gauge fields in compact dimensions with 
(tree-level) brane kinetic terms has been discussed recently in [5| (flat space) 
and [§J [TUJ (warped space) . A particular fermionic case has been addressed 
in jllj . Finally, brane kinetic terms in supersymmetric models have been 
studied in IT3"1 ITU I15j. In all these works, brane kinetic terms are found 
to have a significant impact on model building and phenomenology. And, as 
observed above, they must be generically included for consistency. 

Here we analyse at the classical level the effect of all the possible 0(1/ A) 
brane kinetic terms for scalar, fermion and gauge fields. In particular we 
examine brane terms containing derivatives with respect to the coordinates 
orthogonal to the brane. These terms have not been studied in detail before 1 . 
We find that they give rise to solutions to the equations of motion which ex- 
hibit a rather singular behaviour: they decouple from the brane and become 
insensitive to the magnitude of all brane couplings; furthermore, physical 

1 Except in |16j . where the singular terms were "massaged" via a classical renormaliza- 
tion. We compare this approach with ours in Section [SJ 
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observables for arbitrarily small coefficients of these terms do not approach 
the ones for vanishing coefficients. This is related to the appearance of 5(0) 
singularities at second and higher orders in the 1/A expansion. These singu- 
larities are a common feature of field theoretical models with infinitely thin 
branes [T71 I2]. We also study brane kinetic terms in supersymmetric theories 
and find that supersymmetry requires the presence of higher order kinetic 
terms in the on-shell action. These terms improve but do not cure the singu- 
lar behaviour of the classical fields. More generally, any effective theory will 
contain higher order terms which may significantly change the solutions de- 
rived from the first order Lagrangian. One can then hope that these terms be 
such that the exact solutions are smooth in all the parameters of the theory. 
We will discuss this possibility below. Although most of our considerations 
apply to general field theories in extra dimensions, for definiteness we study 
here theories defined on a flat 5D space with the 5th dimension compatified 
on an S 1 /Z 2 orbifold. Moreover, in order to keep the discussion and calcu- 
lations as simple as possible, we neglect masses and consider kinetic terms 
localized on only one of the fixed points (even if kinetic terms on both branes 
are natural and lead to interesting phenomenology 

The structure of the paper is as follows. In Section|2]we perform the exact 
KK reduction of the free action for particles of spin < 1 in the presence of 
general brane kinetic terms. In Section El we study the supersymmetric case. 
In Section|3]we examine in detail an example elucidating the role of thin brane 
singularities. In Section we argue that classical renormalization makes the 
theory smooth and perturbation theory well defined. This is also interpreted 
as the restriction to a subclass of effective theories with particular reductions 
of couplings. We summarize and discuss our results in Section |B] Finally, we 
collect in Appendix El the general first order results, and give an example of 
renormalization in Appendix [Bl 

2 Exact KK reduction 

In this section we study the exact (i.e., nonperturbative) KK reduction of 
the most general free action with leading order brane terms for massless 
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scalars, gauge bosons and fermions propagating in M 4 x S 1 /Z 2 . We consider 
fields which are functions of the coordinates \x = 0,1,2,3 of M4 and 
y G (— ttR, ttR], of S 1 , and which are eigenstates of Z 2 , i.e., even or odd 
functions of y. For simplicity, we only include brane terms at the fixed point 
y = 0. 

Some of the operators we consider are products of two odd factors times 
a delta function. These terms are sometimes argued to vanish on parity 
grounds. However, these arguments are no longer valid when the (derivatives 
of) wave functions appearing in these terms develop discontinuities across the 
branes [T5j . Deriving and solving the equations of motion when the action 
contains this sort of terms is tricky ^J] > and the same applies to brane terms 
containing derivatives with respect to y. 

In order to treat all these cases rigorously we regularize the branes using 
several (smooth) representations of the delta function, of width ~ e. We 
perform all the calculations with these resolved branes, and only take the 
thin brane limit e — > at the very end. Furthermore, we employ different 
representations in order to check that our results are regularization indepen- 
dent. In some cases we have resorted to numerical methods to solve the 
relevant differential equations. Still, it is possible to see that the numerical 
solutions converge rapidly in the thin brane limit to certain functions, which 
we take as our analytical solutions. The reader can find an explicit analytical 
calculation in the simple case discussed in Section HJ 

2.1 Scalars 

The most general action for a massless complex five- dimensional scalar with 
leading order brane kinetic terms is 2 

S = J d A x J dy {(1 + ad ) dy ~ d y <p ] d y <p 

+ So [\{<P ] d 2 y <P + d^<f>) _ cdytfdy*] } , (2.1) 

2 More generally we could consider a term Sq^^ dy<f) + ^dy(f^(f)), but we take b real 
for simplicity. 
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where 5 = S(y). The dimensionful parameters a, b and c are expected to be 
of order 1/A. In order to perform the KK reduction we expand 

^^ = £#1^)' ( 2 - 2 ) 
„ v 2,tcK 

n=0 

where {f n } is some basis of functions on (—ttR, irR] with the same parity 
character as <fi. Introducing (|2.2j) in the action we find 

S= U*xJ2{ f R (1 + a5 ){^d^d^ 

+ jdy^^^^}, (2.3) 

where we have defined 

= [l + (b + c)5 ]d 2 y + (b + c)5' d y + h -8^ . (2.4) 

In order to diagonalize the quadratic terms, we take {f n } as the set of eigen- 
functions of the following generalised eigenvalue problem: 

Of n = -m 2 n (l + a5o)fn. (2.5) 

The fact that the operator O is hermitian and that the eigenvalues of (|2.5|) 
are nondegenerate, ensures the orthogonality of eigenf unctions with respect 
to the scalar product 

< f>9 >= T~5 1^ ^ i 1 + a6 °) MM ■ ( 2 - 6 ) 
Then, if we also impose, when possible, the normalization 

</„,/„>=l, (2.7) 

the free action reduces to 

S = J d A x {^0 (n)t ^0 (n) - m^Wt^Wj . (2.8) 



The generalised eigenvalue equation (|2.5|) is to be solved with the appropriate 
orbifold constraints. 

When b < or c < 0, it turns out that, for finite (small enough) e, the 
eigenfunctions diverge at some point near y = 0. For this reason, we must 
require b > 0, c > to have a well defined theory. 

Let us now describe the exact solutions in the thin brane limit. For odd 
scalars, the wave functions and KK masses are insensitive to the brane terms. 
The wave functions are 

/ n odd (y) = v / 2sin(m n2/ ), (2.9) 

and the masses, m n = n/R, with n a positive integer 3 . 

For even scalars, the wave functions and masses do not depend on the 
(positive) coefficient c when e — > 0. We distinguish two cases: 

• b = 0. In this case, the solution coincides with the one described in j^J. 
There is a constant massless mode, which depends on the brane terms 
only through its normalization: 

The massive modes are given by 

fn{y) = A n [cos(m n y) - sm(m n \y\)], (2.11) 

where A n is a normalization constant and the masses are the solutions 
of the equation 

tan(7ri?m n ) + = 0- (2.12) 

We observe that if —2itR < a < 0, there is a tachyonic mode, with 
squared mass approaching — 4/a 2 for small |a|/(27ri?). On the other 
hand, the norm squared defined by the scalar product ()2.6|) . ||/|| 2 = 
</,/>, is indefinite. This can be problematic, since zero norm modes 



3 The derivatives of the wave functions, however, have discontinuous values at y = in 
the limit e — > 0, whenever c=^0: f' n (0) = 0=^/^(0 + ). This behaviour at the brane is the 
same as the one for the even component of a fermion, to be described below. 
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have vanishing kinetic term in the reduced 4D action (they cannot be 
normalized as in (|2.7|) ). while modes with negative norm squared have 
the wrong sign in the reduced kinetic term (they are ghosts). We find 
that if a = —2ttR the zero mode has zero norm, and if a < —2ttR it 
has negative norm squared (the massive modes have always positive 
norm). Therefore, we see that for negative values of a the theory is 
either pathological or unstable. 

• b > 0. In this case the solution is independent of all brane terms, 
including b itself (as long as it is strictly positive). We find 

f n (y) = A n sm(m n \y\), (2.13) 

with m n = n+ ft 2 and n integer. In particular, there is no zero mode. 

We see that the KK spectrum and wave functions of even scalars change 
dramatically when an arbitrarily small b is turned on. Note also that the 
wave functions vanish at y = when b > 0. 

2.2 Gauge Bosons 

The Lagrangian for gauge bosons admits brane kinetic terms analogous to 
the ones for scalars, but with the restrictions of gauge invariance. Since we 
are only discussing the quadratic action, we can consider an abelian theory 
without loss of generality. For a gauge boson A M with general 0(1/ A) kinetic 
terms on the first brane, 

S = J d A x J dy {--(1 + aS )F flI/ F^ - i(l + c5 )F 5u F 5 ^ . (2.14) 

The bulk and brane parts respect 5D and 4D Lorentz symmetry, respectively. 
The orbifold projection forces and A 5 to have opposite Z 2 parity. The 
whole Lagrangian is symmetric under the gauge transformations Am — > Am+ 
8m 0, where the gauge parameter 9 has the same parity as A^. Using a gauge 
A 5 = or 0m A M = it is easy to see that the equations of motion are the 
same as the ones for scalars without b terms. At the level of the action, these 
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gauges can be imposed adding the corresponding gauge fixing terms. To 
perform the KK decomposition, however, we choose to work with the gauge 
invariant action. We expand 



oo 



fn(y) 

9yfn(y) A {n), 



4.(*.v) = £^=4°(*), (2-i5) 

n v ZtiK 



n m n V 2irR 

n=0 " 

with = if m = 0, and require 

Of n = -m 2 n (l + a5 )f n , (2.17) 

with 

O = {1 + c5 )d 2 y + c5' d y . (2.18) 

Thanks to the form of the expansion of and the hermiticity of the dif- 
ferential operator O, all the terms get diagonalized. To have canonical 4D 
kinetic terms we also impose the normalization condition ()2.7|) . The eigen- 
functions and eigenmasses can be read directly from the results for (even or 
odd) scalars with 6 = 0. Again, there is no well defined solution for c < 0, 
so we require c > 0. And for a < we have either tachyons or ghosts. 
At the end we are left with the following 4D Lagrangian: 

£ = J2 \-\f^F^ + -d^Al n) d^Af ] + -m 2 n A^A^ + m^d^A^ j . 

n=0 L - - J 

(2.19) 

This Lagrangian is invariant under an infinite set of 4D local transformations: 

4 n) -> A^ + d^ n) , (2.20) 
4 n) ^4 n) +m n 0^, (2.21) 

which is nothing but the KK decomposition of 5D gauge transformations: 

^,y) = f;^eW(x). (2.22) 

n=0 
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The scalars A§ can be interpreted as Goldstone bosons of the spontaneously 
broken gauge transformations 9^ (those with f n (y)=/ constant). They can 
be decoupled from the 4D vector bosons A^ using an gauge. Indeed, 
adding to (|2.19|) the gauge fixing terms 

-^E(^ (n) " + ^ A 5 n) ) 2 , (2-23) 

^ n=0 

we find 

n=0 ^ 

-CmKA^f}. (2.24) 

In the Feynman-'t Hooft gauge, £ = 1, the scalar fields have mass m n , 
whereas in the unitary gauge, £ — > oo, these fields get infinite masses and 
decouple from the spectrum. 



2.3 Fermions 

Let us consider now the most general kinetic action for a five- dimensional 
bulk fermion if) = ipL + ipn, with 7 5 ^L,i? = t4>l,r, and leading order brane 
contributions, 

S = J d 4 x J dy { (1 + a L 5o)^i W>l + (l + a%)^i ^ 

6 b - 

- (1 + -^SojipLdyipR - -5 (d y ip R )ip L 

+ (1 + |5o)MA + ^o(^l)^}. (2-25) 

We have considered again real 6 and c. The two chiralities have opposite 
parity. We take even for definiteness. The KK expansion of the fermions 
reads 

^ L A^y) = ^^=^Lk x )- ( 2 - 26 ) 

10 



In order to diagonalize the 4D kinetic and mass terms, we take the wave 
functions f£ ' R to be the eigenfunctions of the following generalised eigenvalue 
problem: 

OJ* =m n (l + a%)/ n L , (2.27) 

O 2 ti = -m n (l + a R 5 )fK, (2.28) 

where 

Ox =[l + l(b + c)5 ]d y + ^8' , (2.29) 

O 2 =[l+ 1 -(b + c)5o}d y + b -5' . (2.30) 

Because 0\ = —0 2 , the functions f£' R are orthogonal with respect to the 
appropriate scalar products, 

< /, 9 >l,r= ^ dy (l + a L '%) f(y)g(y) . (2.31) 
We also impose the normalization conditions 

< fn, fn >L = < f* , f* >R= 1 , (2.32) 

to render the 4D kinetic terms canonically normalized. 

As in previous cases, we must require b, c > to have well defined solu- 
tions. The solutions then depend on whether b vanishes or not: 

• 6 = 0. There is a chiral zero mode with f R = and a flat left-handed 
wave function, 

fo = f^== ■ (2-33) 

When c = 0, the wave functions for massive modes are 4 

f£ =B n [cos(m n y) + tan(m n 7i\R) sin(ra„|y|)] , (2.34) 
fn =C n [sm(m n y) - a(y) tan(m n 7ri?) cos(m n y)] , (2.35) 

4 Strictly, there is no uniform convergence as e — > due to the behaviour around y = 0. 
We write the limiting functions for \y\ > rj, with rj an arbitrarily small positive number, 
keeping in mind that the interactions of the wave functions with the brane are to be 
calculated before taking the limit e — > 0. 
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with a the sign function and masses given by 



/ oJ' 1 i Tfl 

tan(TrRm n ) + d — tan(Va L a R -y ) = . (2.36) 

These solutions reduce to the ones for even scalar fields when a R = 0. 
Note also that a R is irrelevant when a L = 0. The reason is that only 
when a L forces a discontinuity at the brane is the "odd-odd" a R term 
nonvanishing for e — > 0. When c > 0, the solutions reduce again to the 
ones with a R = 0, i.e., the only effect of a positive c is to cancel the 
effect of a R . 

• b > 0. In this case, the wave functions (for \y\ > r] > 0, see previous 
footnote) and masses reduce to the ones without brane terms: There 
is a chiral zero mode 

fo = 1 , (2-37) 

and a tower of massive modes 

/ n L =B n cos{m n y) , (2.38) 
f R =B n sm{m n y) , (2.39) 

with m n = n/R. In all cases, / n (0) = 0, so that the functions are not 
piecewise continuous. More precisely, 

lim J dy5 e (y)f e L n (y)=0. (2.40) 

Again, we see that the solutions change abruptly when b or c are turned 
on. Moreover, when b > the wave functions vanish at y = 0, so the fermions 
do not couple to brane fields. 



3 Brane kinetic terms in supersymmetric the- 
ories 

In this section we study the impact of brane kinetic terms in supersymmetric 
theories. Besides the intrinsic interest of these theories, one might hope that 
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super symmetry improves the singular behaviour we have observed in the 
previous section. One more reason to look at supersymmetric models is the 
following: we have seen that bosons and fermions with brane kinetic terms 
have different KK wave functions and masses; an intriguing question is then 
how supersymmetry makes the spectra of bosons and fermions in the same 
supermultiplet identical. 

We consider the supersymmetric quadratic actions of massless 5D hyper 
and vector multiplets including all possible 0(1/ A) brane kinetic terms. 



3.1 Hypermultiplet 



Off shell, the 5D hypermultiplet contains two complex scalars, 01,2, a Dirac 
fermion if) (with two Weyl components, ^1,2)5 an d two complex auxiliary 
fields, Fi t 2- We want to construct an action for these fields, invariant under 
5D N=l supersymmetry in the bulk and under 4D N=l supersymmetry 
on the brane. It is easy to write such an action in terms of two 4D N=l 
superfields (201 EE] • The Lagrangian is 



£ — £bulk + £ 



brane 



(3.1) 



with 



and 



A>uik = / d 4 9(^x + $J$ 2 



d 2 e<& x d v <& 2 + h.c. 



£bran e = S S^J d*0 (fll^$l + a 2 ^$ 2 ) 



+ 



j cPBty&xdyQz - C^ 2 d y ^) + h.C. 



(3.2) 



(3.3) 



where, in the y basis, $j = fa + y29ipi + 6 )2 F i . $1 and $ 2 (and their cor- 
responding components) have opposite Z 2 parities. In this case we have 
allowed for complex b and c to keep track of the analytic properties of the 
supersymmetric theory. In components, the bulk and brane Lagrangians read 
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(summation over % is undertood) 

£buik = -<j>p<j>i + id^a^i + F}Fi 

+ [fadyFz + F x d yy2 - fadyfo + h.c.] , (3.4) 

£brane = <5o{o< [ - + id^^^i + F/Fj 

- c[<p 2 d y F 1 + F 2 d y 4 >1 - fadyi/n] + h.c.) } . (3.5) 

The equations of motion for the auxiliary fields are 

[1 + MojFa + [1 + b*6 ]d y 4 + c*d y [5o4] = , (3.6) 
[1 + a 2 6 }F 2 - [1 + c*5 ]d v 4 - b*d v [5 <f>\] = . (3.7) 

We can directly solve them, keeping in mind that the delta functions are 
regularized. Inserting the solutions in the action we find 



C= I -(l + a^o) 



\1 + c5q\ [ 
1 + a 2 5 



-a 



yYl u yYl 



b{ ^ C j ) dy4d y (S 0( f> 1 )+h.C. 

1 + a 2 do 



\b\ : 



1 + a 2 S 

+ 4>Li$ipL + tpRifiipn ~ 4>Ld y lpR + ^ R d y ip L 



(3.8) 



+ S C 



atipLiftipL + a 2 ip R i$ip R + (-bd y ijj R ip L + cip R d y ip L + h.c.) 



We have defined the chiral 4-component spinors ipj, = (^l 0)> ^r = (0 V^)- 
The fermionic sector coincides with the one studied in Section |2] (with b — > 2b, 
c — > 2c). That the spectra of 0i, 02 and ipL,R are identical can be readily 
seen from the equations of motion for the KK wave functions. The equation 
for 0i is 

|1 + (b + c)5 \ 2 d 2 y + [(1 + (b* + c*)5 )(2b + c) + (1 + (b + c)5 )c* 
a 2 |l + (6 + c)5 | 2 lr/Q , L/1 , , „* Nr w, a 2 (5o) 2 



1 + a 2 #o 
-m2(l + ai 5o)(l + a 2( 5o)^ 



]tt + 6(l + (6* + c*)<Jo)(^ 



a 2 <5 



(3.9) 
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and the one for (fr 2 can be obtained from this one by 1 <-> 2 and b <-> c*. For 
the fermions, 

{ [1 + (6* + c*)5 ]^ + c^}/* = m n (l + ai 6 )tf, (3.10) 
{ [1 + (b + c)5 ]d y + b5' }ti = -m n (l + a 2 8 )f* (3.11) 

It is then trivial to show that iterating the fermionic equations we obtain the 
scalar ones and thus the spectra of fields in the same 4D N=l supermultiplet 
are exactly the same, as implied by supersymmetry. We see that in order for 
the scalar solutions to reproduce the fermionic ones, higher order (singular) 
terms are required in the scalar sector. They are reminiscent of the famous 
5(y) 2 terms of 

This example illustrates the fact that higher order operators can modify 
radically the spectrum, although supersymmetry does not prevent the sin- 
gular behaviour we observed in the last section for the coefficients b and c. 
Nevertheless, we can consider the theory with 6 = and c = 0, which be- 
haves smoothly and looks stable under quantum corrections, as these terms 
are protected by nonrenormalization theorems. Note, however, that higher 
order D-terms, which can in principle be induced, may reintroduce the sin- 
gular behaviour. 

3.2 Vector multiplet 

The 5D off-shell vector multiplet consists of a 5D vector Am, two Weyl 
gauginos Ai^, a real scalar E, a real auxiliary field D and a complex one F x . 
Under D=4 N=l supersymmetry they form a vector supermultiplet, V, and 
a chiral one, x EDI EE] : 

V = -Oa^OA^ - i^OXi + idWXx + h 2 9 2 D , (3. 12) 

X= -^(Z + iA 5 ) + V2i9\ 2 + 6 2 F x , (3.13) 
v2 

where V is written in the Wess-Zumino gauge and x i n the "y" basis. Since 
we are interested in the free action we consider an Abelian theory. The gauge 
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transformations are 



V -> V + A + A f , 
X -> X + >/29 y A , 



(3.14) 
(3.15) 



with A a chiral superfield. The gauge invariant Lagrangian including 0(1/ A) 
arbitrary brane terms is C = £tmik + £branc with 

£ bu ik = -\J d 2 ew a w a + h.c. + J d 4 e (d y v - -±=( x + x f )) , (3.16) 



and 

^branc Oq 



a I d 2 9W a W a + h.c. + c / d 4 6 ( d y V - -j=(x + X f ) 



(3.17) 

The corresponding action is 5D N=l supersymmetric in the bulk and 4D 
N=l supersymmetric on the brane. V and x have opposite Z 2 parity. In 
components the Lagrangian reads 

£buik = - l -F MN F MN + l -d^d^+ l -D 2 -W y D + F x Fl 
- Xiia^d^Xi - X 2 icr^d fl X 2 + A 2 <9j / Ai + A 2 <9j / Ai , 

^brane 



(3.18) 

- a -F^F, u - °-F^F„ + |^E^E + a -D 2 - cZd y D + cF x F\ 



aXxia^d^Xi - cX 2 ia^d^X 2 + cX 2 d y X 1 + cX 2 d y Xi] . 



(3.19) 



Using the equations of motion for the auxiliary fields, F x = and (1 + 
aS )D = —dy[(l + c5o)Y]i we find the following Lagrangian for the dynamical 
fields: 

C = ~(l + a8 )F^F^ + ^(l + c8 )F^F tl5 

- Ul + c5 )d^d^ - l-J—d y [(l + c5 )nd y [(l + c5 )S] (3.20) 

Z Z 1 + OOg 

+ (1 + a5o)X L iftX L + (1 + c8 o )\ R i0\ R + (1 + c5 )(X R d y X L + d y X L X R ) , 

where we have defined the chiral 4-component spinors A^ = (Ai, 0) and A^ = 
(0,A 2 ). 
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Again, the equations of motion for bosons and fermions in the same 4D 
supermultiplet are identical, which leads to the expected degeneracy of the 
spectra. Indeed, the relevant differential equations are 



8y (1 + C5 )dy /, 



■A 



(3.21) 



n 



for the gauge boson, 



d 2 y (1 + cW* 




dy (l+C^)/ n 5 



m 2 n (l + a5 )fl 



(3.22) 



for the scalar, and 



dyfn = -mnl 



R 



m n (l + a5 )f„ : 



(3.23) 
(3.24) 



for the fermions. If we now iterate the first order fermionic differential equa- 
tions we find that they are identical to the ones of their bosonic counterparts. 

Using the results in the previous section we observe that the KK masses of 
all the fields are the same as the ones for gauge bosons. In particular, they do 
not exhibit a Mike behaviour 5 . This can be understood as the combination 
of gauge invariance, N=l supersymmetry relations and the fact that both 
chiralities combine to produce 4D Dirac masses. On the other hand, there is 
no finite solution when c < 0, and the wave functions of even E and Xr are 
forced to vanish at y = when c > 0. 

To summarize this section, supersymmetry improves the behaviour of the 
KK masses and wave functions, but it does not cure it. 



4 Singular solutions, the thin brane limit and 
perturbation theory 



We have seen that the brane operators which contain derivatives with respect 
to y give rise to solutions behaving in a quite singular way. In this section we 
5 For even this is so in the fermionic sector because supersymmetry fixes a R = c. 
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study in greater detail a simple example which exhibits the essential features 
of the general cases studied above. 

Consider the following free Lagrangian of a 5D massless fermion: 

£ = 4) L i$ipL + i*Ri $ipR - i>LdyipR + ip R dyil) L 

-\WLd y ^R + {d y i> R )i) L ]. (4.1) 

We take ipi even and ipR odd. The coefficient of the brane term, b, is naturally 
of order 1/A. Once more, we assume that the delta function is regularized. 
The equations of motion can then be derived in the usual way. They read 

{l + ^5 )d y ^ R = i^ L , (4.2) 
[{l+ b -6 )d y + b -5' }^ L = -z^ R . (4.3) 

These equations admit a zero mode solution: ipL,R(x,y) = f L ' R (y)£,L,R{x), 
with $£,l,r = 0. This solution is chiral since the boundary conditions for the 
odd component, 4>r{Q) = ^R^R) = 0, together with ()4.2|) with a vanishing 
rhs, imply that / = everywhere. This leaves a homogeneous first order 
equation for f L with no further constraint. Using the following representation 
of the delta function, 

we readily find 

= 2Ntt- y —^- (4.6) 

be + 27c{y 2 + e 2 ) V ' 

with iV a normalization constant. The first equation, ()4.5|) . holds for any 
regularization of the delta function. Observe that when b < and e < \b\/2ir 
in ()4.6j) . the function is singular at the points y = ±a/— e(6/2 + 7re)/7r; more- 
over, its integral is divergent. Therefore we have to restrict b to nonnegative 
values. The limit e — > then depends only on whether b = or b > 0. When 
b vanishes the limit of (J4.5j) is simply f L (y) = N. For b > 0, f L approaches 
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N everywhere except around y = 0, where it goes to zero in such a way that 



does not couple to fields living on the brane. Moreover, it is insensitive to 
other possible brane terms and to the precise (positive) value of b itself. 

Obviously, in the thin brane limit, the zero mode is not smooth in the 
coefficient b: the solution for a small positive value of b differs drastically 
(around the brane) from the one for vanishing b, and there is no finite solution 
for any negative b. On the other hand, for a thick brane with e > |6|/27r, the 
zero mode is well defined and smooth in b. 

This behaviour also shows up in the limit of the perturbative expansion 
in 1/A of the finite e result (J4.5j) : 



We see that the breaking of perturbation theory is caused by the functions 
8q, n > 1, which are singular when e — > 0. Strictly, the expansion (J4.7|) 
is only valid for small \b\/e, and therefore it is incompatible with the limit 
e — > 0. However, the situation is analogous to the one in standard (renormal- 
izable or nonrenormalizable) quantum field theories, where the quantum UV 
divergences make the loop expansion badly defined. Nevertheless, in these 
theories it is well known that one can work order by order after renormal- 
ization: each order is finite and (at weak coupling) smaller than the lower 
ones. Since our classical divergences have essentially the same origin as the 
usual quantum ones (the unknown UV of the theory), it is natural to apply 
the same renormalization procedure to obtain well defined, finite predictions 
at each order in perturbation theory. We develop this point of view in the 
next section. 

5 Renormalization and critical theories 

The divergences we have encountered in perturbation theory are actually a 
common feature of field theories with infinitely thin defects. They appear, 



J dy5(y)f L (y) — ► 0. Therefore, for positive b the zero mode of the fermion 




(4.7) 
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for example, in classical electrodynamics with point sources. These singu- 
larities indicate a breakdown of the effective theory at scales in which the 
details of the fundamental theory — for instance, in our case, a finite thick- 
ness of the brane or stringy effects — cannot be neglected. In this sense, they 
are analogous to the usual UV divergences in radiative corrections. In |2*2*] . 
Goldberger and Wise have shown how to deal with classical divergences in 
brane theories of codimension greater than one by the usual renormalization 
procedure of quantum field theory. The same idea can be applied to brane 
theories of codimension one with more singular operators: adding appropri- 
ate counterterms it is possible to cancel the 5(0) singularities order by order. 
Let us see how this can be achieved to 0(1/ A 2 ) in the example of the previ- 
ous section. In this simple case it is sufficient to add to the Lagrangian the 
following counterterms: 

C ct D -j6 2 J L i0ip L - j5 2 {$ L d y ij) R + h.c.) . (5.1) 

Then, it is straigthforward to check that, to second order, the zero mode 
solution is well defined: 

/ L = iV(l-^ + O(l/A 3 )), (5.2) 

and the kinetic term in the action is finite: 

f R dy (1 - ^8 2 ) (f L ) 2 = N 2 (2nR - b + 0(1/ A 3 )) . (5.3) 

This example is quite trivial, but it illustrates the main properties of the 
counterterms: 

1. They are operators of higher dimension. This is related to the fact that 
the brane coefficients are dimensionful and the theory is nonrenormal- 
izable. 

2. The divergences are not pure numbers, but functions of y localized at 
the branes, i.e., functions like 8$. Notice that also one loop corrections 
in an orbifold without brane terms require the introduction of coun- 
terterms in a nontrivial background: the delta functions at the fixed 
points. 
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3. The coefficients of these singular operators are determined (in a given 
regularization) by the requirement that the divergences must be can- 
celled. This is analogous to the fact that the infinite parts of usual 
counterterms are fixed by the structure of the radiative corrections. 
However, in our case there is no freedom to adjust at will the finite 
part of the coefficients, as the operator itself is "infinite" . This has a 
very important consequence: 

4. The coefficients of the higher dimensional operators with more than 
one delta are not independent parameters of the effective theory. They 
do not require new observables to be fixed by experiment. 

Of course, a full calculation at a given order will include not only the clas- 
sical part we have considered so far, but also loop diagrams and, possibly, 
contributions from higher order operators in the effective action. All these 
contributions may include new thin brane singularities, which should be can- 
celled. This can be achieved adjusting the coefficients of the counterterms 
accordingly. In particular, they will contain quantum infinities. In this way, 
the usual quantum divergences and the thin brane singularities are subtracted 
all at once. 

It should be noted that the renormalization procedure we propose is not 
identical to the one in [22] . where the classical divergences are just numbers 
which can be absorbed into the bare first order couplings. In fact, the 5D 
theories we are analising differ from the 6D model in [22j * n a ^ least two im- 
portant aspects. First, we find power-law, instead of logarithmic divergences. 
Second, we are dealing with operators with a richer structure, due to the in- 
clusion of derivatives orthogonal to the branes. This seems to require in the 
counterterms nontrivial backgrounds with more than one delta function. 

We would also like to comment on the similarities and differences of our 
approach to the one in JH]. There, the singularities in 5D theories with 
brane operators, like the ones we are discussing, were analised in perturba- 
tion theory. The authors proposed to substitute the singular brane terms 
by terms without y derivatives, giving the same solutions to the equations of 
motion outside the core of the brane. The resulting action was interpreted as 
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a renormalized action. However, we observe that this "massage" (we follow 
the terminology of JE]) by itself does not lead to finite solutions since in the 
original action the divergences near the branes induce divergences also for the 
solutions far from the branes, and even for observables such as KK masses. 
In the "massaged" theory this divergent behaviour comes from explicit diver- 
gences in the coefficients of the nonsingular operators, which translate into 
divergencies in the boundary conditions just outside the core of the branes. 
On the other hand, subtracting the explicit divergences in the mentioned 
coefficients leads to really finite results. This can be achieved by adding 
counterterms such as the ones in (j5.1j) . which has the additional advantage 
that the behaviour inside the brane is also well defined. We should remark 
that, in their application of "massaging" to the study of the stability of the 
compact Randall-Sundrum model, the authors of [15] considered arbitrary 
finite coefficients of operators without y derivatives. Hence, they did use a 
renormalized theory. 

In Appendix B we work out the renormalization of a more involved ex- 
ample: the propagator of a scalar field in the presence of a "6" brane term. 
A more complete study of renormalization of theories with brane terms will 
be presented elsewhere. 

This renormalization process can be interpreted in a complementary way. 
We have seen that, in general, theories with brane kinetic terms have a sin- 
gular behaviour in the thin brane limit (before renormalization), such that 
perturbation theory is spoiled. However, there is a subclass of "critical" 
theories which are well behaved and predictive: the ones with a tower of 
operators identical to the counterterms introduced above (plus all the opera- 
tors present initially). The coefficients of these operators are not independent 
parameters: both their finite part and the part depending on the quantum 
renormalization scale are functions of the coefficients of bulk operators and 
brane operators with only one delta. In this sense, these critical effective 
theories correspond to a particular reduction of couplings [23 121] • In fact, 
such relations should exist in the low energy effective theory if the funda- 
mental theory is to be well behaved for infinitely thin branes. For instance, 
in [T7] it was hoped that M-theory would directly give well defined quantities 
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without the 6(0) that were observed in the supergravity description of the 
Hofava-Witten model. 6 A particularly simple example of critical theory (at 
least at the classical level), in which all the coefficients of operators of order 
higher than one vanish, is the one for a gauge boson with c = 0. Indeed, 
in this case one can directly perform the KK reduction to all orders in a, as 
done in [3] and in Section |21 

6 Conclusions 

We have studied 5D theories for spins < 1 with general brane kinetic terms. 
We have shown that some terms which had been neglected in the past play 
an active role, especially in the case of fermions. We have first calculated 
the KK masses and wave functions in these theories to all orders in the coef- 
ficients of the brane terms, using a combination of analytical and numerical 
methods. We have found that some terms change smoothly the solutions 
without brane terms, while others — the ones with derivatives with respect 
to y — change them abruptly, destroying the perturbative hierarchy of the 
effective theory. We have also built supersymmetric free actions with brane 
terms for the hyper and vector multiplets and studied their KK decomposi- 
tions. Supersymmetry does not solve in general the difficulties produced by 
the second kind of operators, although it can alleviate them in some cases. 
We have then discussed in detail a particular case, showing that the singular 
behaviour is caused by the presence of infinitely thin branes. If we insist on 
performing perturbative calculations, we have to face <5(0)-like divergences at 
0(1/ A 2 ) and higher orders. At this point, at least three different paths can 
be followed. 

First, we can simply take the exact results at their faith value and ac- 

6 It might be objected that supersymmetric theories cannot be critical — or, equivalently, 
renormalization breaks supersymmetry — since in our discussion above supersymmetry de- 
termined completely the higher order operators of the on-shell theory, and we have seen 
that they do not eradicate the singular behaviour. This is actually not so, since we in- 
cluded only first order operators in the off-shell action. Higher order off-shell operators 
compatible with supersymmetry can give rise to a critical theory. Of course, the on-shell 
supersymmetry transformations will be different from the ones without the new operators. 
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cept a theory with such behaviour. From a phenomenological point of view, 
the possible realistic models would be constrained by the decoupling of the 
wave functions from the branes and, in the nonsupersymmetric scenario, the 
absence of scalar zero modes. In particular, bulk fermions would not cou- 
ple to a boundary Higgs, while a nonsupersymmetric bulk Higgs would not 
have the required zero mode. Even more problematic is the fact that per- 
turbation theory breaks down. It seems very difficult, if not impossible, to 
construct a predictive interacting field theory in extra dimensions without 
some hierarchy controlling the size of the operators. 

Second, we can work with thick branes instead of zero width branes 
(see for recent calculations with kinetic terms on thick branes). This 
is natural if the brane is a domain wall. The disadvantage of this approach 
is that the results depend on the substructure of the brane and are thus very 
model dependent. Furthermore, we would like to be able to define sensible 
field theories in orbifold compatifications with vanishing brane width. 

Third, we can renormalize the theory to render it free from both classi- 
cal and quantum divergences in perturbation theory. This is equivalent to 
working with a "critical" theory containing a tower of operators whose co- 
efficients satisfy certain relations. Here we have just argued that this is a 
sensible approach to construct extra dimensional models with thin branes. 
We have not described the renormalization program in detail nor proven that 
it can be carried out to all orders. This subject is under investigation. 

Including mass terms in the bulk and on the branes may be relevant for 
phenomenology, but it does not alter in an essential way the results given 
here. The same holds for brane kinetic branes on the second brane. Finally, 
in a complete model one should also take into account brane interactions, 
which may require the introduction of new counterterms. 
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A First order solutions 

In this appendix we present the KK expansions to first order in 1/A. The first 
order results are directly finite, but they are only meaningful when the theory 
is renormalized, so that higher orders are also finite. Perturbatively, the 
solutions to the generalized eigenvalue problem can be computed in different 
ways. If the full analytic solutions are known for a smooth regularization 
of the brane, we can expand them in powers of the dimensionful coefficients 
a, b, c and then take the limit e — > 0. Alternatively, we can perform the 
perturbative expansion from the very beginning in the eigenvalue equations 
and normalization conditions and solve them order by order (again for a finite 
e and taking the limit of thin brane only at the end of the calculation). 

We find that the first order KK masses are invariant under a <-> —b for 
scalars and under a L <-> —b and a R <-> — c for fermions, while the wave func- 
tions are invariant up to possible delta functions. This can be understood in 
the following way. In the scalar sector, one can perform the field redefinition 

*-(l-^o)6 (A.l) 

which to first order only has the effect of redefining the coefficients 6^0 
and a — > a — b. In the fermion sector, the field redefinitions 

- (l-^i, (A.2) 
■4>r - (l-^o)fe (A.3) 

redefine b — > 0, a L — > a L — b and c — > 0, a R — * a R — c. All these field 
redefinitions generate higher order operators as well. These are responsible 
for the singularities observed at higher orders. 



25 



Scalars and Gauge Bosons 

Consider an even scalar with general values of a, b, c (gauge bosons satisfy 
the same equations with the condition 6 = 0). The solution to first order is, 
for the zero mode, 

A <"> = 1 + l=£-5*' (A.4) 
and for the massive modes, 

f n (y) =V2cos(m^y)-^=So 

+ ^ b ^[cos(m^y) + 2m^(7rR-\y\)sm(m^\y\)}, (A.5) 

with rrin^ = n/R, and masses 

Note that the "odd-odd" parameter c has no effect to this order. 
For an odd scalar there is no zero mode and the massive ones are 

f n (y) = V2sm(m^y) 

+ [sinrfy) - 2m®a{y){*R ~ \v\) ^s(m^y)] , (A.7) 

with a(y) the sign function and masses 

m n = -(l + -^—) . (A.8) 

The "odd-odd" parameters, which are now a and b do not contribute to this 
order while the "even-even" parameter c has a non-trivial effect, in contrast 
with the results in Section EJ 



Fermions 

For the even component (which we take here to be the LH component), this 
case is identical to the scalar case except for the term with a R , which, just 



26 



like c, does not contribute to first order. The wave functions and masses to 
0(1/ A) are 

f L (y) = 1 + - 5*0, (A.9) 

for the zero mode (with no zero mode for the odd chirality), and 

f^y) = V2cos(m^y)-^=5 

+ ^= b -^[cos(m^y) + 2m^(nR- \y\) sm(m^\y\)], (A.10) 

fK(y) = V2sm(mWy) (A.ll) 

+ ^f^|-[ sin ( m n 0) V) -2m^a(y)(7cR- M)cos(mf y)], 

(A.12) 

m " = fl( 1 + wj- (A13) 

for the massive modes. 



B Renormalization of the scalar propagator 

Here we give another example of classical renormalization. We consider the 
scalar theory (j2.1j) with a = c = for an even scalar, and calculate per- 
turbatively the classical propagator to 0(1/ A 2 ), i.e., 0(b 2 ). We work in 
momentum space for the first four coordinates and in position space for the 
5th coordinate. The brane delta functions £1X6 , cLS usual, assumed to be reg- 
ularized and this time we indicate this explicitly through a subindex e, to 
differentiate these deltas from the ones which do not need regularization. In 
perturbation theory, the propagator can be written 

A(p; yi, m) = A (p; Vl , y 2 ) + Ax(p; y u y 2 ) + A 2 (p; y h y 2 ) + 0(b 3 ) , (B.l) 

with A n proportional to b n . The zeroth order propagator satisfies the equa- 
tion | 

[d 2 yi + p 2 } A (p; Vl , y 2 ) = -- [5( Vl - y 2 ) + 5{y x + y 2 )} , (B.2) 
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with boundary conditions 

[d yi A (p; Vl , ife)] w=0 + = , V2^0 , (B.3) 

[^AoCpjyi,^)]^^- =0, yrf-KR. (B.4) 

The solution is 

a / \ cos(py<) cos(p(ttR - y> )) 

Ao(p 1/1,^2) = ^—7 — ^ , (B.5) 

with y< = Min{|j/i|, 1 2/2 1 } , y> = Max{|j/i|, 1 2/2 1 } - The first order propagator is 
given by the zeroth order propagator with one insertion of the brane kinetic 
term. It reads 

b f *~ 
Ai(p;yi,jfe) = g / dz 5 e(z)A (p;y 1 ,z)(d 2 +d 2 )A (p;z,y 2 ) 

2 P 2 A (p; 0)A (>; 0,2/2) - <J(yi)A (p; yi, y 2 ) 
-Aofa^jfcWs/a)]. (B.6) 



In the second line we have used (IB.2j) and taken the limit e — > 0, which is 
well defined. The singularities arise at the next order. The second order 
propagator has two insertions of the brane kinetic term: 

b 2 f *~ 
A 2 (p;2/i,y 2 ) = - d^idz 2 6 e (z 1 )8 e (z 2 )Ao(p;y 1 ,z 1 ){ d\ +d Zi )A (p;z 1 ,z 2 ) 

x(S 2 Z2+ d 2 Z2 )A (p;z 2 ,y 2 ). (B.7) 

Using (|B.2|) for the internal propagator and integrating on z 2 , we find a 
contribution 

b 2 f r 

[A 2 (p;?/i,2/ 2 )] div =j / dzi |p 2 5 e 2 (2;i)Ao(p;2/i,zi)Ao(p;2;i,?/ 2 ) 

+ <9 Z1 [5 e (z 1 )A (p;y 1 , z x )]d zi [5 e (z 1 )A (p; z u y 2 )] 

- 5 2 { Zl ) [Ao(p; 2/1, *i) ( +<) A (p; y 2 ); 

(B.8) 
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which diverges when e — > 0. To renormalize the theory to second order (at 
the classical level) we add the following counterterms to the Lagrangian: 

Cct = T { ~ 5 e(y) d ^ 1d ^ - 9 v HvW] d y [Uy)<f\ + <E(v) {^d 2 y <p + <9^V) } . 

(B.9) 

Then, the divergences are cancelled and taking e^Owe find the finite result 

Ar(p;y 1: y 2 ) = j [2p 2 A (p; y u 0)+5( yi )] A (p; 0, 0) [2p 2 A (p; 0, y 2 )+5(y 2 )] . 

(B.10) 

We have checked that the counterterms (jB.9|) also render the KK reduction 
finite to order 0(1/ A 2 ). The explicit delta functions in the propagator can 
be translated through a field redefinition into brane interaction terms. For 
points yi and y 2 away from the first brane, the renormalized classical prop- 
agator reads 

A rcn (p; Vi, 2/2) =A (p; 2/1, 2/2) - bp 2 A (p; Vl , 0)A (p; 0, y 2 ) 

+ &yA (p; 2/1, 0) A (p; 0, 0) A (p; 0, y 2 ) + 0(b 3 ) 

[cos(p y< ) + |psin(py<)] cos (pjirR - y>)) + ^ 3 
2p[sm(p%R) — |pcos(p7ri?)] 

(B.ll) 

Actually, working out the renormalization to all orders and resumming the 
series, it can be shown that the exact renormalized propagator is equal to 
the explicit expression in the last line of (jB.ll|) 26J. The KK masses and 
wave functions can be easily obtained from (jB.ll|) : they are the poles and 
the square roots of the residua of the propagator. They agree with the ones 
for an even scalar when b = c = and a = —b. In fact, the explicit expression 
in (jB.ll|) is identical to the exact propagator of the theory with b = c = 
and a = —b. As we have observed in Appendix El the field redefinition (jA.l|) 
eliminates the b operator from the action (j2.1J) and redefines a — ► a— 6. It also 
introduces new operators at higher orders, which give rise to singularities. 
The result in (jB.ll|) shows that performing the field redefinition ()A.1|) in the 
action with the counterterms (jB.9|) . the second order terms cancel out (at 
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least when a = c = initially). Conversely, performing the field redefinition 
and neglecting the higher order operators which are generated amounts to a 
renormalization of the theory. The delta functions in (jA.l|) account for the 
ones in A rcn . 
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